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Abstract 



1 Introduction 



Sigma models arise frequently in particle physics and condensed-matter physics as low- 
energy effective theories. In this paper I compute the exact free energy at any tempera- 
ture in two hierarchies of integrable sigma models in two dimensions. These theories, the 
SU(N)/SO{N) models and the 0(2P)JO(P) x 0(P) models, are asymptotically free and 
■^j- ■ exhibit charge fractionalization. When the instanton coupling 9 = n, they flow to the 

CO \ SU(N)i and 0{2P)\ conformal field theories, respectively. I also generalize the free en- 

ergy computation to massive and massless perturbations of the coset conformal field theories 
g ! SU{N) k /SO{N) 2k and 0{2P) k /0(P) k x 0{P) k . 

o 

^— > , 

Two-dimensional sigma models have been the subject of a huge amount of study because they are 
interesting toy models for gauge theories, because they often arise in experimentally-realizable 
condensed-matter systems, because this is the highest dimension in which they are naively renor- 
malizable, and because of the powerful theoretical methods applicable. 

One of the nice things about sigma models is that the same model can often describe com- 
pletely different physics. The reason is that in many situations, the precise sigma model of interest 
follows mainly (or sometimes entirely) from the symmetries. For example, sigma models often 
arise in theories of interacting fermions invariant under some group G. If some fermion bilinear 
gets an expectation value manifestly invariant under some subgroup H, then the excitations at 
low energy can be described by a field taking values in G/H. Put another way, the expectation 
value gives the fermions mass at some scale M. One can then integrate out fermionic excitations, 
leaving only bosonic G/H excitations with masses below M. The sigma model describes the 
interactions of these low-energy excitations, and is independent of many of the details of the 
original theory. This is why vastly different theories may end up having the same low-energy 
physics. 

Two-dimensional G/H sigma models all have a global symmetry group G, even though the 
fields take values in the smaller space G/H. This is one big difference between two and higher 
dimensions. In higher dimensions, the symmetry G of these sigma models would be spontaneously 
broken to H, and in the effective low-energy-theory, the G symmetry is not manifest. In other 
words, in higher dimensions the sigma model describes the physics of the massless Goldstone 
bosons. However, the Mermin-Wagner-Coleman theorem says that in two dimensions continuous 
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symmetries cannot be spontaneously broken. The way these sigma models satisfy this theorem 
is to give the would-be Goldstone bosons a mass and keep the original global symmetry intact. 

In particular, many interesting sigma models in two dimensions are asymptotically free. At 
large energies the interactions are weak, but at low energies the interactions are strong. Naively, 
there seems to be no mass scale in the theory (the coupling constant g is dimensionless), but a scale 
/j, appears in the theory as a result of short-distance effects which need to be renormalized. The 
coupling g depends on this scale. At \x large, <?(//) is small, so the theory is effectively free, while 
as [i decreases, g(n) increases. In renormalization-group language, there is an unstable trivial 
fixed point at g = 0. For G/H sigma models, the manifold G/H has dimension dimG— dinrff, so 
as g — > the theory reduces to dimG— dimH free bosons. 

Very elaborate techniques of perturbation theory have been developed to describe sigma 
models in the regime where g(fi) is small (see ||]). However, when a sigma model is being 
used as an effective theory, it is only applicable to the relevant physics at low energies, where 
-C M. Usually in this regime, g{y) is large. Thus while the perturbative techniques give 
valuable information, they may not tell the whole story. To understand the regime where g(fj,) 
is large, one must utilize alternative techniques. Large- N expansions are a common and useful 
tool. However, for most applications N is small. For example, an application of great current 
interest in the condensed matter community is in sigma models describing disordered systems. 
These sigma models are derived by using the replica trick, which requires sending N — > at the 
end of the computation. Obviously, large-A expansions are not necessarily going to be reliable 
here. 

Luckily, for two spacetime dimensions there are other non-perturbative methods applicable. 
Many sigma models are integrable, with an infinite number of conserved currents. The resulting 
conserved charges constrain the system, making exact computations possible, even at strong 
coupling. The aim of these paper is to attempt to discuss a number of aspects of integrable sigma 
models. I will derive the exact free energy at finite temperature and in the presence of a magnetic 
field. This makes it possible to compute the susceptibility and specific heat. It also makes it 
possible to understand exactly the effects of the theta term, a modification of the sigma model 
action which drastically changes the low-energy physics. 

One extremely interesting question is if g(fM) continues to increase as fi decreases, or if it 
reaches a fixed point. The existence of a fixed point obviously affects the physics enormously. In 
the sigma models describing disordered systems, g is related to the conductance of the system. 
If there is a fixed point, the system is a conductor, with conductance determined by the value 
of g at the fixed point. If there is no fixed point, the system is an insulator. In the former case, 
the excitations of the model are massive, while in the latter, they are massless. For the models 
discussed in this paper, a non-trivial fixed point appears if a theta term is added to the sigma 
model action. The theta term has no effect on perturbation theory. Nevertheless, as shown in 
|§, H, H, its presence can result in the appearance of a fixed point at large g, completely unseen 
in perturbation theory. 

There are two sets of sigma models to be discussed in this paper. Their actions can be written 
conveniently in terms of a symmetric matrix field <& as 

S = -tr J d 2 x (1) 

along with the constraint 

$t$ = = j (2) 
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where I is the identity matrix. The constraint (]2|) that <I> be unitary can easily be imposed by 
adding a potential like Atr (<^<3? — I) 2 and taking A large. In theories with interacting fermions, 
this often results from introducing a bosonic field to replace four-fermion interaction terms with 
Yukawa terms (interactions between a boson and two fermions). Integrating out the fermions 
then gives such a potential for the bosons and hence the sigma model. 

In the first set of models discussed in this paper, the field takes values on the SU(N)/SO(N) 
manifold. This corresponds to taking $ to be a symmetric, unitary N x N matrix of determinant 
1. The simplest case, N = 2, corresponds to the manifold SU{2)/ SO{2) being a two-sphere. This 
is because a general symmetric unitary 2x2 matrix of determinant one can be written as 

V 1V3 v\ - iv 2 J 

where v%, v 2 and Vs are real and obey (t>i) 2 + (V2) 2 + (V3) 2 = 1- 

In the second set of models discussed in this paper, the field takes values on the 0(2P) /O(P) x 
0{P) manifold. This corresponds to taking $ to be a symmetric, orthogonal, real, and traceless 
2P x 2P matrix. There are several correspondences between the two sets of models, because 
50(6) = 5*7(4)/Z 2 , 50(4) = 517(2) x 5C/(2)/Z 2 , and 50(3) = 5C/(2)/Z 2 . The case P = 2 
therefore reduces to two decoupled copies of the two-sphere, whereas the sigma model with 
P = 3 is equivalent to the 5f/(4)/50(4) sigma model. 

The reason these G/H manifolds can be described in terms of symmetric matrices is as follows. 
In both cases, the global symmetry G acts on the field <3? as 

$ -» U$U T (3) 

where U is a unitary matrix of determinant one. This transformation preserves the fact that # is 
a symmetric matrix with determinant ±1. In the 0(2P)/0(P) x O(P) sigma models, the matrix 
<1> is also real. To preserve this reality, U must be real as well, so G = 0(2P). The eigenvalues of 
a orthogonal matrix must be ±1, and if the matrix is traceless as well, there must be the same 
number of +1 and —1 eigenvalues. The field $ in this case can diagonalized with an orthogonal 
matrix U, so $ can be written 

$ = UAU T $ G 0(2P)/0(P) x O(P), 

where U is in 0(2P), and A is the matrix with P values +1 and P values —1 on the diagonal. 
Different U can result in the same <3?: the subgroup leaving $ invariant is H = 0(P) x 0(P). 
This is why the space of symmetric orthogonal traceless matrices is indeed 0(2P)/0(P) x O(P). 
For the SU(N)/SO(N) models, U can be any unitary matrix of determinant one, so the global 
symmetry G is indeed SU(N). Field configurations here can be written in the form 

§ = UU T $ G SU{N)/SO{N) 

where U is in SU(N). The subgroup H leaving $ invariant is SO(N). For example, $ = I for 
any real U in SU{N), i.e. if U is in the real subgroup SO(N) of SU(N). This is why H = SO(N) 
here. 

Under renormalization, the matrix $ preserves its form: e.g., it remains symmetric. In other 
words, the space G/H preserves its "shape" under renormalization, with only the overall volume 
changing. The effect of renormalization is to increase the curvature (increase g). These sigma 
models are all asymptotically free, so going to high energies decreases g. This behavior happens 
for all sigma models on symmetric spaces G/H (where H is a maximal subgroup of G). 
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With the action (|T|), there is no fixed point at large g. However, if one adds a theta term, there 
is a non-trivial fixed point in these sigma models ||] . A theta term affects field configurations with 
non-zero winding number n, which are called instantons. The winding number is a topological 
invariant; roughly speaking, it counts the number of times the field configuration wraps around 
the two-dimensional spacetime. The theta term is then 

S e = inO. (4) 

If the winding number n takes integer values, the theory is periodic under shifts of the coupling 
9 to 9 + 2tt. This is why the coupling 9 is often called an angle. However for the general cases 
considered here, n can take just two values, and 1. This means that 9 takes just two values 
here: 9 = and 9 = tt. The variables n and 9 should be thought of as Fourier conjugates. Adding 
the 9 term to the action amounts to doing a discrete Fourier transform. 

For the sphere sigma model (the case N = 2 or P = 2 here), n takes integer values. It 
was argued in j|, ^] and proven in [|j that when 9 = tt in the sphere sigma model, there is a 
non-trivial fixed point at large g. This behavior is widely believed to persist in other models with 
a 9 angle (see |?J for a review) . An important question is therefore whether the existence of these 
non-perturbative fixed points in sigma models at 9 = tt can be generalized. In Q, it was shown 
that the SU(N)/SO(N) and 0(2P)/0(P) x O(P) sigma models have non-trivial fixed points 
at 9 = tt. The former fixed points are described by the SU(N)\ WZW theory, while the latter 
are described by the 0(2P)\ WZW theory. The exact spectrum and S matrices were found, and 
used to compute the energy at zero temperature in the presence of a background field. This 
computation essentially proves the existence of these fixed points. 

It is the purpose of this paper to complete this proof by studying the behavior of these models 
at finite temperature. I will compute a c-function ||] which clearly shows how the field theory 
flows from the trivial fixed point {g = 0) to the non-trivial fixed point at some large value of g. 
This computation also makes it possible to compute the specific heat and susceptibility at both 
9 = and 9 = tt, a fact which will be useful in other work |J. 

In section 2, I discuss the thermodynamic Bethe ansatz formalism necessary to do the com- 
putation. In section 3, I compute the free energy at any temperature for the massive 9 = sigma 
models. In section 4, I compute the free energy for the massless 9 = tt models. In section 5, I 
discuss some related coset models. I conclude in section 6 by discussing the symmetries of these 
sigma models, and the prospects for generalizing these results to other sigma models. 



2 The Thermodynamic Bethe Ansatz 

The proof that the sphere sigma model has a non-trivial fixed point at 9 = tt utilizes the in- 
tegrability of the model at 9 = and tt |||, 1C]. Integrability means that there are an infinite 



number of conserved currents which allow one to find exactly the spectrum of quasiparticles and 
their scattering matrix in the corresponding 1 + 1 dimensional field theory. The quasiparticles 



for 9 = are gapped and form a triplet under the SU{2) symmetry [11], while for 9 = tt they 
are gap less, and form SU(2) doublets (left- and right-moving) This is a beautiful example 
of charge fractionalization: the fields (^1,^25^3) form a triplet under the SU(2) symmetry, but 
when 9 = tt the excitations of the system are doublets. To prove that this is the correct particle 
spectrum, first one computes a scattering matrix for these particles which is consistent with all 
the symmetries of the theory. From the exact S matrix, the c function can be computed. It was 
found that at high energy c indeed is 2 as it should be at the trivial fixed point, while c = 1 as 
it should be at the SU(2)i low-energy fixed point Q. 
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As an even more detailed check, the free energy at zero temperature in the presence of a 
magnetic field was computed for both 8 = [1^] and ir JO]]. The results can be expanded in a 
series around the trivial fixed point. One can identify the ordinary perturbative contributions 
to this series, and finds that they are the same for 8 = and tt, even though the particles 
and S matrices are completely different This is as it must be: instantons and the 8 term 
are a boundary effect and hence cannot be seen in ordinary perturbation theory. One can also 
identify the non-perturbative contributions to these series, and see that they differ. Far away 
from the trivial fixed point, the non-perturbative contributions dominate and cause a non-trivial 
fixed point to appear when 8 = tt. The computation of the energy at zero temperature in a 
background field was done for the SU(N)/SO(N) and 0(2P)/0(P) x O(P) sigma models in §. 

In this paper I will compute the exact free energy at any temperature, and thus compute the 
c function. I will use a technique called the thermodynamic Bethe ansatz (TBA), which I will 
describe in this section. 

2.1 The exact S matrix 

An integrable field theory possesses an infinite number of conserved currents and charges. The 
symmetries strongly constrain the dynamics, but without making the system trivial. The con- 
straints are why the theory is "solvable". In this context, solvable means that some quantities 
can be computed exactly. These constraints imply that once the particle spectrum is known, the 
exact S matrix can be found. Integrable models have the striking property that in a collision all 
momenta are conserved individually, and that the n-body S matrix factorizes into a product of 
two-body ones. This two-body S matrix is completely elastic, meaning that the momenta and 
energy of the particles are conserved individually, not just overall. Internal quantum numbers 
can change in a collision, so the S matrix is not necessarily diagonal. There are two possible ways 
of factorizing the three-particle amplitude into two-particle ones; the requirement that they give 
the same answer is the Yang-Baxter equation. There have been hundreds of papers discussing 
how to solve this equation, so I will not review this here. For a detailed discussion relevant to the 
sigma models here, see e.g. 

0, HI 0i El- Solutions arising in the sigma models will be given 

below. 

One of the useful characteristics of having particles in representations of a Lie algebra is 
that their S matrix can be written in terms of projectors onto representations of this algebra. 
The invariance of the G/H sigma model under the Lie-group symmetry G requires that the 
S matrices commute with all group elements. The S matrix can then be conveniently written 
in terms of projection operators. A projection operator V a maps the tensor product of two 
representations onto an irreducible representation labelled by a. By definition, these operators 
satisfy VaPb = $abPb- Requiring invariance under G means that the S matrix for a particle in 
the representation a with one in a representation b means that the S matrix is of the form 

s ab m = ^2ffmv c (5) 

c 

where /3 = j3 a — is the difference of the rapidities, and the f£ b are as of yet unknown functions. 
The sum on the right-hand side is over all representations c which appear in the tensor product 
of a and b; of course ^2 C V C = 1- In an integrable theory, the functions fc b (P) are determined by 
requiring that the the S matrix satisfy the Yang-Baxter equation. 

I define the prefactor F ab (f3) to be the coefficient f£ b in (|5|) where the highest weight of the 
representation c is the sum of the highest weights of the representations a and b. The Yang-Baxter 
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equation does not give this prefactor. To obtain it, one needs to require that the S matrix be 
unitary, and that it obey crossing symmetry. With the standard assumption that the amplitude 
is real for (5 imaginary, the unitarity relation S^(f3)S(f3) = I implies S((3)S(—f3) = I. The latter 
is more useful because it is a functional relation which can be continued throughout the complex 
(3 plane. Crossing symmetry is familiar from field theory, where rotating Feynman diagrams by 
90° relates scattering of particles dj and bj to the scattering of the antiparticle aj with bj. 

Multiplying any S matrix by function F((3) which satisfies F((3)F{—(3) = 1 and F(in — 0) = 
F{(5) will give an S matrix still obeying the Yang-Baxter equation, crossing and unitarity (this is 
called the CDD ambiguity). To determine F((3) uniquely, one ultimately needs to verify that the 
S matrix is consistent with the bound-state structure, and that it gives the correct c function. 



2.2 Fusion 

In this paper, I derive the TBA equations for the sigma models by utilizing fusion. Fusion is a 



method of finding new solutions of the Yang-Baxter equation from known ones [16|. One starts 
with a solution where the states are in some representation of a symmetry algebra. Then one can 
find new solutions in other representations, just as one takes tensor products of representations. 
The usual place fusion appears in the study of exact S matrices is in what is called the bootstrap 



(see e.g. [14]). In many integrable models, various particles can be thought of as bound states 
of other particles. The bootstrap procedure relates the S matrices of bound state to those of its 
constituents. However, fusion is a more general procedure than just the bootstrap. It can be 
used to relate S matrices of different models. This fact will prove very useful here, because when 
the S matrices are related, the TBA equations are related as well. This observation enables the 
computation of the TBA equation for integrable sigma models. 

Formally speaking, fusion relies on the observation that at certain values of 0, the coefficients 
of some of the projectors in the S matrix vanishes. This means that some particles can be treated 
as being composites: they are composed of "constituent" particles at specific rapidities. I avoid 
calling the composite particles bound states, because this implies that the composites and the 
constituents are both particle states in the same theory. This is the not case in general. For 
example, the only particles in the sine-Gordon model at 2 = 8ir particles are in the spin-1/2 
representation of SU(2), while in the sphere sigma model, the only particles are in the spin-1 
representation of SU(2). Fusion means that the S matrices are related, even though the theories 
are different: the spin-1 particles are composites of the spin-1/2 ones. 

I will demonstrate fusion in theories with SU(N) symmetry. The two-particle S matrix for 
two particles in the A^-dimensional vector representations of SU(N) contains two terms: one 
involving the projector Vs onto the symmetric representation, the other Va onto the antisym- 
metric representation. This is because the tensor product of two symmetric representations in 
SU(N) decomposes into the irreducible symmetric (N(N + l)/2 dimensional) and antisymmetric 
representations (N(N — l)/2 dimensional): 

(N) ® (AT) = (N(N - l)/2) (N(N + l)/2). 

For SU(2), the antisymmetric representation is the singlet, so this statement means that two spin 
1/2 representations tensored together is the sum of the spin-0 and the spin-1 representations. 
The vector- vector S matrix for SU(N) is determined by requiring that it satisfy the Yang-Baxter 
equation. It is 

S VV m=F VV (P)(v S + l±^V A ). (6) 
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The function F vv (/?) is the prefactor I defined above. It must be consistent with unitarity, 
crossing and the bootstrap. A "minimal" solution of these constraints means the S matrix has 
no poles in the region < Im{(3) < tt. The minimal solution here is 

jpWfaS _ r ( 1 ~^) r (^ + ^) ( , 

r min\l J ) / a \ / a 1 \ V'J 

r (l + A rf-^ + i 



For a given model, the prefactor F ab ((3) may or may not be the minimal solution. This prefactor 
is crucial to the physics, but the fusion procedure is valid for any F ab {(3). 

At (3 = —2iri/N, S in (|6|) involves only the projector onto the symmetric representation. 
The fusion procedure means that particles of rapidity (5s in the symmetric representation can be 
treated as being composed of two constituents in the vector representation, of rapidities (3s—iTr/N 
and Ps + iw/N. The reason this works is described in [16]. The variable j3 in the S matrix is 



the difference of the rapidities of the two particles, so when (3 = 2m /N, the antisymmetric 
combination is effectively projected out. The Yang-Baxter equation ensures that this projection 
survives any scattering. In other words, if two vector particles are in the symmetric combination, 
they can scatter from other particles and change state. However, if their rapidity difference is 
2m/N, the final state of these two particles will still be part of the symmetric representation. 

Because particles in the symmetric representation are composed of vector constituents, the 
S matrices are related as well. The S matrix for scattering two particles in the symmetric 



representation has three terms. In the language of weights [17|, the symmetric representation has 
highest weight 2/ii, and the tensor product is 



(2/ii) <g> (2/ii) = (4//i) (2 W + n 2 ) © (2/i 2 ) 



The S matrix is 



qSS (fi) _ F SS (fi) ( v i P+MK y + 2m/N (3 + 4vri/iV \ 

S W-F {(3) + - _ 4m/N + p _ 2m/N p _ Am/N T^m ) ■ (8) 

The explicit form of the projection operators is given in p|. The minimal solution of the unitarity 
and crossing constraints F^ n {f3) has no poles in the region < Imf3 < tt, and is 

r(i-£)T(J J + 2± 



,SS ra\ I 3 - 27Ti / N V ~ V^" 17 



Note that F%f n (f3) differs from F^ n ((3 + 2m/N)(F^J n ((3)) 2 F^Y n ((3 - 2m /N); the prefactor does 
not automatically follow from the fusion procedure. 

In cases where the composites are bound states of the constituents (all are particles in 
the same theory), then the bootstrap procedure relates the prefactors of composite scatter- 
ing to those of constituent scattering. However, the fusion does not make such a require- 
ment in general: the prefactor F ((3) does not necessarily follow from F (/3). All the fu- 
sion procedure does is determine the overall form of the S matrix for the composite parti- 
cles and ensure that it obeys the Yang-Baxter equation. Although one might expect that 
F SS (P) = F vv ((3 + 2iTT/N)(F vv {(3)) 2 F vv ((3 - 2m /N), I will show that below this is not true 
in general here. In another words, the CDD ambiguity may be resolved in different ways in the 
constituent and composite theories. 
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2.3 The free energy of an integrable theory 

Once the exact S matrix is known, the exact free energy as a function of mass, temperature, and 



magnetic field can be computed by using the thermodynamic Bethe ansatz (TBA) [18, |19|| . This 
enables one, for example, to compute thermodynamic quantities like the susceptibility. It also 
allows a very substantial check on any assumption of integrability. The reason is that at a critical 
point, the free energy is known exactly - it is related to the central charge of the corresponding 



conformal field theory [20]. Thus the free energy computed from the TBA must give this result 
in the limit where the mass of the particles goes to zero, and the system is at the unstable UV 
fixed point. 

The TBA requires a relation between the density of states of the particles to the actual particle 
density. This relation is called the Bethe equation. If the particles are free, this is trivial: the 
density of states is independent of the particle density. If the scattering is completely elastic and 
diagonal, this relation is easy to derive. This is because a diagonal two-particle S matrix is the 
boundary condition the phase shift in the wave function: 

ip(xi,x 2 ) = e *w*i-HM*a forxi«x 2 

ip{x 1 ,x 2 ) =e ipiX1+ip2X2 S( Pl ,p 2 ) forzi»x 2 (10) 

In a state of J\f particles, the Bethe equation follows by requiring that one-dimensional space of 
length L be periodic, and that the wavefunction be invariant under sending any of the coordinates 
X{ — > Xi + L. First consider the case where there is only one kind of particle in the spectrum, 
with two-particle S matrix S((3\ —fa). The requirement of periodicity of the wavefunction 
ip(xi, X 2, ■ ■ ■ xj\f) yields the relations 

AT 



eimsinhe.LJJsiPi- (3j) = 1 (11) 



One can think of this intuitively as bringing the particle around the world through the other 
particles; one obtains a product of two-particle 5-matrix elements because the scattering is 
factorizable. This is the generalization of the free-particle momentum quantization condition 
p = 2nir/L. 

The Bethe equation is written in terms of the density of states P{(3) and the density of 
rapidities p((3). The former is defined so that the number of allowed states with rapidities 
between (3 and (3 + d(3 is P(/3)d[3, while the number of states actually occupied in this interval 
is p(/3)d(3. The quantization condition relates the two. Taking the derivative of the log of (|Tl| ) 
yields 

/oo 
df3'<S>(f3- P')p(f3') (12) 
-oo 

where $(/?) = j-^s InS(P). This is easily generalized to the situation where there is more than one 
particle in the spectrum, as long as the scattering is diagonal. Let S a b be the S matrix element 
for scattering a particle of type a from one of type b. Defining densities P a and p a for each type 
of particle, the Bethe equations are 

/oo 
d(3' $ ab (P - (3')p b (/?')• (13) 

where 



S 



Once the Bethe equations are known, the TBA equations and hence the free energy can be 
derived. This is done by minimizing the free energy, using ( |l3| ) as a constraint. The result is 
most conveniently written in terms of the "dressed particle energies" e a (f3), defined by 

EM = I (14) 

For simplicity, I have set all chemical potentials and background fields to be zero. The resulting 



TBA equations are [0, [19] 

e a {(3) =m a cosh/?-^T f° f^<M/? - /?') In (l + e~^' T ) (15) 

For free particles, & a b = and the e a just reduce to the particle energies. The form of the TBA 
equations reflects the fact that in all integrable particle theories of this type, it is either proven 
or assumed that the particles fill levels like fermions: at most one particle in a level. The free 
energy per unit length F is given in terms of these dressed energies e a . It is 

F(m, T) = —T m a J ^cosh/?ln(l + e- £a(/3)/T ) (16) 

In the IR limit m a —* oo, the gas of particles becomes dilute, and interactions can be neglected. 
The free energy becomes 

lim F(m,T) = -TV m a f°° — cosh cosh(/3)/T (17) 

a J ~°° 

This integral can be done, yielding a Bessel function. 

Calculating the free energy using the TBA allows an extremely non-trivial check on the exact 
S matrix. In the limit of all masses going to zero, the theorem of |^] says that the free energy 
per unit length must behave as 

vrT 2 

lim F = —c uv (18) 

m a -*0 D 

where cjjv is the central charge of the conformal field theory describing this UV limit. The 
number cjjv can usually be calculated analytically from the TBA, because in this limit the free 



energy can be expressed as a sum of dilogarithms |21]. The cuv computed from the TBA must 
of course match the cjjv from the field theory. This provides an extremely non-trivial check not 
only of the S matrix, but of whether the entire spectrum is known. All particles contribute to 
the free energy, so if some piece of the spectrum is missing or if an incorrect particle is included, 
the correct cuv will n °t be obtained. 

The TBA computation is much trickier if the scattering between particles is non-diagonal, as 
is the situation for the models of interest here. The Bethe equation is much harder to derive, 
because as one particle is going around the periodic world, it can change state as it scatters 
though the other particles. This requires introducing the "transfer matrix" T for bringing the 
a given particle through the others; since the scattering is not diagonal, the final state is not 
necessarily the same as the initial. To define T explicitly, I first introduce the scattering matrix 
T a b((3) for bringing a particle of type a and rapidity through M particles and ending up with 
a particle of type b. Thus the different T a b make up a set of s 2 x matrices, where s is the 
number of different types of particles. The scattering is completely elastic, so the rapidities do not 
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change even though the scattering is not diagonal. This means T a b(P) depends on the rapidities 
f3\ . . . (5j\f as well as f3. Let S a b^ c d(Pi — (h) be the two-particle S matrix element for scattering 
an initial state a(/3i)6(/?2) and ending with a final state of c{j32)d{(3\). Then the components of 
Tab can be written in terms of the S matrix elements as 

(T ab (m . . . Mtlt:*% = E s ^d lfl ip - 0i)s flC2 -> d2h 09 - /%) . . . s fMCN ^ dNh {p - p M ) 

where the sum is over the intermediate states f\ = 1 . . . s, f 2 = 1 . • • s, . . . , fj\f = 1 . . . s. The 
matrix T follows by exploiting the fact that all the S matrices of interest at zero relative rapidity 
just permute the colliding particles. In other words, 5o&-»cd(0) = —$ac&bd- Thus setting (3 = (3 a 
effectively turns the oft 1 particle so that it scatters through all the others. This is precisely what 
is needed for the TBA. To put periodic boundary conditions on the system, one sums Taa over 
all a. The result is that 

T((3 Q \(3 U ...fa)=Y, Taa(P = Pa\Pa+l,- • • PM, Pi, • • • , Pa-l)- (19) 
a 

This is a s-^^ 1 x s-^^ 1 matrix. 

The TBA requires finding the eigenvalues A(P a \Pi, ■ ■ ■ Pm) of T. The crucial effect of the S 
matrix satisfying the Yang-Baxter relation is that the T((3) commute for different f3. This ensures 
that T{(3) can be simultaneously diagonalized for all f3 by a /3-independent set of eigenvectors; 
only the eigenvalues depend on (3. The quantization condition (11) is generalized to 

e imaSinhA * L A(/3 a |/?i, ...p Af ) = l (20) 

This must hold for all particles a = 1 . . . N . In the limit of large M, A depends on the particle 
densities instead of the individual rapidities. Henceforth I will just write A(/3). For the cases of 
interest here, finding the eigenvalues A(/3) is quite difficult, but has been done in 22, ^3], |24|| . 
The Bethe equations are still of the form (|l^) , and the TBA equations are still of the form ( |l5| ) . 
However, extra particles, known as "pseudoparticles" or "magnons", enter the equations. These 
particles appear in the equations just as if they were a particle species, but with m a = 0. I will 
give examples of the explicit form of these equations below. 

The transfer matrix has very nice properties under fusion, because the fused S matrices are 
products of the constituent S matrices. The case of most interest here is when particles in the 
representation with highest weight [i a are fused to give particles in the representation 2fi a . Then 
the transfer matrices for A^/2 fused particles is related to the product of transfer matrices for J\f 
constituents. The reason it is a product is because both constituents must be brought around 
the world in the fused transfer matrix. The precise relation is 

T 2 ^(p a \/3i,...l3 M/2 ) = C(p a )T^(p a + r J \l3 1 + r l ,f3 1 -r},...l3 M/ 2 + V,pAr/2-v) x 

(A, - n\(3i + rj, Pt - v, ■ ■ ■ Pn/2 + V, Pn/2 ~ rj) (21) 

The rapidity difference of the constituents is 2rj. The reason for the extra factor C(/3) is that the 
prefactors of the S matrices need not satisfy the exact fusion relation, as discussed above. This 
constant of proportionality is 

C(3) = TT F 2 »°"°(/3-/3 a ) 

a=i F^ aVa (/9 — Pa + rj)F^ al/a (/3 — /3 a — 7]) 

where the particle with rapidity p a is in the representation v a . Given this relation between 
transfer matrices, the eigenvalues obey the relation 

A 2 ^ (p) = C(P)A fla (p + rj) A^ (/3-T7). (22) 
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3 Massive sigma models 



In this section I will derive the TBA equations for a variety of massive sigma models. I start 
with the sphere sigma model, before going on to the more complicated cases. 

3.1 The sphere sigma model 

One of the best-known sigma models is the sphere sigma model, where the field takes values 
on a two-sphere. In the G/H language I have been using, this corresponds to G = SU{2) or 
G = SO(3), and H = U(l) or H = SO(2). The TBA equations were derived originally by taking 
the limit of certain integrable fermion models [po] , p6| , and conjectured on different grounds in 
[27]. I will rederive the TBA equations here directly from the S matrix, because this is the 



method which generalizes most simply to the more general sigma models of interest. 

In a two-dimensional G/H sigma model, the global symmetry group is G. Therefore the 
symmetry group of the sphere sigma model is G = SO(3): the symmetry corresponds to rotations 
of the sphere. The particles of this model were shown long ago to be in the spin-1 representation 
of SO (3) Their S matrix was derived by solving the Yang-Baxter equation directly, and is 
given by (8) with N = 2 and 

Since this S matrix is non-diagonal, one needs to diagonalize the transfer matrix as described 
in the last section. The way to do this is to first solve the problem for particles in the spin-1/2 
representation of SU(2), and then use fusion to find the answer for spin 1. For particles in the 
spin-1/2 representation of SU(2), the two-particle S matrix is given by (0) with N = 2. This S 
matrix is four-by-four, since there are just two different kinds of particles (spin up and down). 
The choice 

pVV _ r?VV 
r N=2 — r min 

gives the S matrix of the sine-Gordon model at the coupling /3 2 = 87r in the usual conventions. At 
this coupling, the dimension of the cos j3(f> perturbation is two, so that it is marginally relevant; 
the £7(1) symmetry of the sine-Gordon model is enhanced to SU(2). Another name for this model 
is the SU(2) Gross-Neveu model. 

For particles in the spin-1/2 representation of SU(2), the Bethe equations were derived 70 
years ago, in the original paper by Bethe himself 2Hj. The reason is that the transfer matrix for 
the spin-1/2 representation of SU{2) as defined in (l{j) precisely corresponds to the transfer matrix 
of the Heisenberg spin chain. In the limit of large number of particles jV, the eigenvalues of the 
transfer matrix follow by adopting the "string hypothesis". This means that the eigenvalues A(f3) 
of the transfer matrix defined in ( |l9| ) are expressed in terms of densities pk((3), with k = 1 . . . oo. 
These are the pseudoparticles discussed above: they enter the TBA equations as if they were real 
particles with no mass term. (I have somewhat abused the conventional notation: most authors 
would not use the ~ here, but it makes subsequent relations less confusing.) The other density 
entering the equations is the density of particles po(/3). This is the total particle density, with 
contributions of both spin up and spin down particles. 

Bethe's result for the eigenvalues is 

, oo 

- In AO?) = * po(/3 ) + £ ff («0 * fttf) (24) 
ap 3=1 
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where convolution integrals are denned as 



/*<?(/?)= / dp'fip-p'W). 



The kernels are given explicitly in the Appendix. The kernel Y( N > comes from the prefactor of the 
S matrix. This only affects the coupling to the total particle density, and not the pseudoparticles, 
because it contributes an overall factor Y\^ = iP VV (P — Pa) to the transfer matrix. Now I can 
write down the first of the Bethe equations, by taking the derivative of the log of (p0|). This gives 



oo 

2vrP (/3) = mcosh/3 + * Po {P) - * pj(p). (25) 

j=i 

where m is the mass of the particles. Pq is the total density of states for the particles. The other 
Bethe equations relate the densities of states for the pseudoparticles to particle and pseudoparticle 
densities. They are 

oo 

27T Pj (/3) = * po(J3) - £ 4° } * p0) (26) 



1=1 



where the density of string states Pj is 



p j = Pj + Pj 



Note that all the Bethe equations are of the form (13), with no mass term for the pseudoparticles. 

Using identities in the appendix, all the Bethe equations (including that for Pq) can be written 
in the compact form 

2nPM = 5 j0 m cosh + J^' —^pz^pi^') ( 2? ) 



Here the indices j and I in the incidence matrix Ij^ = Sjj+i + run from 0,1,..., oo. 

Note that the right-hand-side involves the hole densities, not the particle densities. This Bethe 
equation is conveniently represented by the diagram in figure 1. With these equations, it follows 
from the standard TBA calculation that the TBA equations ( p^|Jl^ ) hold, with 



and 



r(°°) 

cosh(/3) 



rrij = 5jom cosh j3. 

These equations were first derived in the context of the sine-Gordon model at P 2 — ► 87r in |2^] . 
One can easily check that the free energy has the correct properties. In the UV limit m/T — ► 0, 
one obtains the correct central charge cjjv = 1 from (|l^). This follows from a now-standard 



analysis, involving expressing the free energy as a sum of dilogarithms (see e.g. [21, 19, ^]). In 
the IR limit, the generalization of ([17]) to the case with pseudoparticles is 

F = mT(l + e- fl H) 1/2 f° ^ cosh/3e — cosh(/3)/T 
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Figure 1: The incidence diagram for the SU (2) Gross-Neveu model (the sine- Gordon at j3 2 -> 8tt). 
The circles represent the functions e a ; the filled node represents the fact that the equation for eo 
has a mass term. The line represents the coupling between the functions in the TBA equations. 



For particles with rrij ^ 0, e~ e ^ ca " > vanishes. However, the pseudoparticles have no mass term, 
and here one finds that e - ^ 00 ) = (j + l) 2 — 1 for j > 1. This means that the free energy in the 
IR limit is that of 2 types of particles of mass m, as it must be. 

It is now simple to get the S matrices and TBA for the sphere sigma model by using fusion. 
The fusion procedure says that the spin-1 particles in the sphere sigma model can be viewed as 
having the spin-1/2 particles as constituents. As explained above, a spin-1 particle (in a represen- 
tation with highest weight 2 P i) is composed of a pair spin-1/2 particles (each in a representation 
with highest weight P i) with rapidities Pi + in/ 2 and Pi — in/2. The transfer matrix for ftf /2 
spin-1 particles is related to that for the M spin-1/2 particles by the relation (|2l|) with r\ = in/2. 
Because the two transfer matrices are related in this way, the Bethe equations for the sphere 
sigma model follow from those above after a few modifications. The eigenvalue of the sphere 
sigma model transfer matrix follows from the spin-1/2 eigenvalue (|24|), and the fusion equation 
(H). It is 

, oo 
jL lnA sphere {p) = Z ^ * p (P) + * Pj (P) (28) 

where 

r^\p) = af°\p + in/N) + o { °° ] {p - in/N) 

with N = 2 here. The first term in ( p8| ) arises from the prefactor of the sphere S matrix (|23"|), 
with 

-O-"**. - FT?- 

The explicit expressions for and are given in ( |70|) and ( |67|) in the appendix. Using this 
expression for the eigenvalue in ( |20| ) gives 

oo 

2vrP 09) = m cosh p + Z<® * Po (P) - £ * Pj (/?) . (29) 

i=i 

The Bethe equations for the densities of states of the pseudoparticles (^) are modified because 
the real particles come in pairs with rapidities P ± in/2. Thus for the sphere sigma model 

oo 

2n Pj (P) = * po(P) - 4^ * PiO 3 ) ( 3 °) 

1=1 

for j > 1. 

By using the identities in the appendix, the Bethe equations ([?(], ^) can be put in the unified 
form 



2nP J {P) = ^omcosh^ + ^X)^ / ^7^)^ 

1=0 J ~°° 



(31) 
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Figure 2: The incidence diagram for the sphere sigma model 



The indices j and I here run from . . . oo. Above, the incidence matrix 1^ was associated with 
SU(s). Here, the incidence matrix 1^ is associated with 0(2s): iff = 25 ji — C^ 2s \ where 
C°( 2s ) is the Cartan matrix for 0(2s). Explicitly, 

Zji = s j,i+i + &3,i-i + 5j,2$i,o + 5jfl$i,2 ~ Sj,i8i,o ~ 8jfl8i,i (32) 

This Bethe equation is conveniently represented by the diagram in figure 2. 

With these equations, it follows from the standard TBA calculation that the TBA equations 
(ffMD hold, with 



j(oo) 

cosh(p) 

and 

rrij = 5jom cosh f3. 



One can easily check that the free energy has the correct properties |27|. In the UV limit 
m/T — > 0, one obtains the correct central charge cjjv = 2 by the standard dilogarithm analysis. 
In the IR limit, one finds that 



F = mT(l + e - e2 ^ V2 f°° ^ cosh 



-mcosh(/3)/T 



As with the spin-1/2 system, the functions obey e <e j( 00 ) = (j + l) 2 — 1 for j > 1. This means 
that the free energy in the IR limit is that of 3 types of particles of mass m, the spin-1 triplet. 

3.2 SU{N) Gross-Neveu models 

To find the sigma model free energy, it is best to first perform the analysis for the vector particles 
and then use fusion. The appropriate field theory with particles in the vector representation of 
SU (N) is the SU(N) Gross-Neveu model (also sometimes called the chiral Gross-Neveu model) 



|3l| , 32]. Its similarities and differences with the sigma model were discussed at length in ||. The 
vector particles in the SU(N) Gross-Neveu model have the S matrix (||). The prefactor Fq^(P) 
is not the minimal one given in ([?]). It is instead 

where 

_ sinh(i(/? + 2^/AQ) 
X{(3) ~ sinh(I(/?-2™/A0)" (33) 
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Note that X = 1 for N = 2, so the sine-Gordon model at f5 2 — > 87r is indeed the SU(2) Gross- 
Neveu model. 

The pole at /3 = 2iri/N in this factor X((3) means that for N > 2, the vector particles in the 
Gross-Neveu model have bound states in the antisymmetric representation. Upon completing 
the bootstrap procedure, one finds that the model has bound states in all the antisymmetric 
representations with a indices, a = 1 . . . N — 1. These are called the fundamental representations, 
and they have highest weight p a . The particles can be expressed as bound states of a particles 
in the vector representation. These have mass 

m n = m sinh ( — — 



N 

The representation with highest weight pN~a is the conjugate of the representation p, a , because 
of the invariant e tensor. For example, the N representation has highest weight /xjv-1 anci mass 
m-Tv-i = mi. The bootstrap procedure gives the S matrices for all scattering of these particles. 
The scattering is not diagonal, but it is diagonal in the representation labels. When a particle 
in representation a with rapidity (5 a scatters, the final particle with rapidity f3 a must be in some 
state in same representation a. This means that the two-particle S matrix prefactors can be 
labelled by F ab . The vector- vector prefactor F vv = F 11 in this new notation. The explicit 
prefactor Fq N is necessary for the calculation, and is given in (|6tf ) in the appendix. 

Computing the Bethe equations for the SU(N) Gross-Neveu models looks extremely difficult 
or impossible. Remarkably, the computation has already been done in ^4|] by using fusion. 
Here the Bethe equations are found for any simply-laced Lie algebra G, when the particles are 
in any representations with highest weight m/ij where pi is a fundamental weight of G, and m is 



an integer. This work was generalized to non-simply-laced groups in 33]. The fusion procedure 
gives functional relations like (|2l| ) for all the T a ((3 a \f3\, . . . ,f3j\f) |23; ]- The label a here indicates 
that the a particle is in the representation with highest weight p a . These functional relations 
relate various T a . The prefactors F ab ((3) need to be computed, but the explicit S matrix is not 
needed: all the relevant physics is contained in the representation theory and in the fusion. From 
the functional relations and a few mild analyticity assumptions, the eigenvalues of T a and the 
Bethe equations can be derived in the limit of a large number of particles. 

The Bethe equations for the general case require the introduction of pseudoparticle densities 



and densities of states into the Bethe equation (13). Here the pseudoparticle densities p a j and 
densities of states P a ,j(P) are labelled by two indices. (In the literature, this is usually called 
a nested Bethe ansatz.) The index a runs from 1 to N — 1 for SU(N). For the N = 2 case 
treated above, this index takes only one value can be suppressed. The index j is the same index 
as before, running from 1, ... oo for the pseudoparticles. The functions p a ,o and P a fl are defined 
respectively as the density and density of states for all the particles in the representation \i a . 
It is consistent to define separate densities for each representation, because the particles cannot 
change representation when scattering. For all values of a and j, P a j = p a j + p a ,j- 

The computation of the TBA equations directly from the SU (N) Gross-Neveu model S matrix 
was done in [ 34 1 . The eigenvalues of the transfer matrix T a are |23|, ^4| 



N-l oo 

In A a GN (f3) = £ Y<p * Pb> oW) + E ^ * P«AP) (34) 

P 6=1 3=1 

where the kernels are given explicitly in the Appendix. The kernel comes from the prefactor 
F ab of the S matrix. It couples the density of states of real particles in representation a to the 
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density of particles in representation b. The first of the Bethe equations follows from ( |20| ) , and is 

JV-1 oo 

2vrP Q , (/?) = m a cosh (3 + £ Y$ } * p 6>0 (/?) - £ of * p oJ (/?) . (35) 

6=1 j=l 

The other Bethe equations relate the densities of states for the pseudoparticles to particle and 
pseudoparticle densities. They follow from [23, as well, and are 

N-l oo 
6=1 Z=l 

where P a j = p a j + p a j . Explicit expressions for these kernels are given in the Appendix. Note 
how all these equations reduce to those in the last subsection by setting N = 2. 

By using the fact that A and K are inverses, and the identities in the appendix, all the Bethe 
equations (35,36) can be written in the combined form [34] 

N-l 00 

2npaj (/?) = S j0 m a cosh - £ E * A% * Ph>l {fi) (37) 

6=1 1=0 

Here the indices j and / run from 0, 1, . . . , 00. With these densities, the dressed energies e a j(/3) 
are defined as in fll4|). It follows from the standard TBA calculation that the TBA equations 



(gjyEg) hold, with 
and 

m a j = 5jom a cosh (3. 

The TBA equations can be rewritten in a much more elegant form by using the fact that A and 
K are inverses, and the simple relation between K and the incidence matrix 

ljf° = 5j, ,_j + <^+i j, Z = 1 . . . A - 1 (38) 

The result is 

^ " ^ a6 7-oc 2vr2cosh(iV(/3-/?0/2) in l 1 + e J 

^ A lnfl + e-^C ')) (39) 



00 „ c 

In, •/ — < 



J_ O0 2*2caBh.(N{p-P')/2) 

This is a substantial simplification because the equation for e a j only involves "adjacent" functions 
e a ,j±i and e a ±ij. These equations are displayed schematically in figure 3. The dashed and 
unbroken lines account for the different minus signs in (|39"|). Note that the masses do not appear 
in rewritten TBA equations (|39|), although they appear in the original ones. When using the 
form (|39|), the asymptotic conditions 



e a ,o(P — > 00) — > m a cosh/3. 

must be imposed. 

This free energy of the SU(N) Gross-Neveu model has the correct properties. In the UV limit 
m/T — ► 0, one obtains the correct central charge cj/y = A — 1 from the dilogarithm analysis. In 
the IR limit, one finds that each representation contributes one term to the free energy, with the 
correct multiplicity (e.g. N for the vector representation a = 1, N(N — l)/2 for the antisymmetric 
representation a = 2). 
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Figure 3: The incidence diagram for the SU(N) Gross-Neveu model. There are N — 1 rows and 
an infinite number of columns. 

3.3 SU(N)/SO(N) sigma models 

Here I find the TBA equations for the SU(N)/ SO(N) sigma model, generalizing the analysis for 
the sphere sigma model, which corresponds to N = 2. The TBA analysis is related to that for 
SU(N) Gross-Neveu models via fusion for all N. 

The SU(N) I SO(N) sigma models have a Lagrangian description (||) in terms of a a symmetric 
and unitary matrix field. The particles of the sigma model are in all representations with highest 
weight 2fi a , a = 1 ... iV — 1 ||. The representation with highest weight 2/j,\ is the symmetric 
representation. The two-particle S with both particles in the symmetric representation is given 
by (H|) with prefactor || 

F SS (P) = X(f3)Fif n {(3) 



where the minimal factor is given in (^), and X(j3) is in (|33]). The pole in X(j3) at /3 = 2m/N 
means that particles in the representation 2^2 are the bound state of two particles in the symmet- 
ric representation 2/xi. Because the factor X{f3) is the same as that of the SU(N) Gross-Neveu 
model, the masses are the same: 

m a = msui(ira/N) 

for the sigma model as well. However, the multiplicites are different because the former are in 
representations with highest weight /x a , while in the latter they are in representations with highest 
weight 2fj, a . 

As discussed above, f3 = —2-ki/N, the S matrix (B) is entirely in the symmetric channel. 
Therefore, the particles in the symmetric representation 2/xi can be viewed as composites of those 
in the vector representation n\. The same is true for all the particles in the SU(N)/ 'SO(N) sigma 
model: those in the representation 2\x a are composites of two particles in the \i a representation. 



Because of this relation between S matrices, the transfer matrices are also related by (21) 23]. 
This means that the resulting TBA systems are closely related, and all the densities are labelled in 
the same way. Explicitly, the Bethe equations for the SU(N)/SO(N) sigma model are obtained 
from those of the Gross-Neveu model by two modifications. The kernel coming from the S 
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matrix prefactor is replaced with while the kernel crj^ is replaced with Tj\ defined by 

rj s) ((3) = af ((3 + m/N) + af ((3 - m/N) . (40) 



The sigma model version of (^) is 



N-l oo 



2vrP ai0 (/3) = m a cosh /? + ^ * p b , (/3) - ]T rj^ * p a>j (p). (41) 

6=1 j=l 

while the Bethe equations for the pseudoparticles are 

N-l oo 

2irp ad ((3) = * Pa , (f3) -EE A t ] * #2° * WW ( 42 ) 

6=1 {=1 

Explicit expressions for these kernels are given in the appendix. Note how all these equations 
reduce to those of the sphere sigma model by setting N = 2. 

The different kernels in the Bethe equations of course mean that the TBA system is not quite 
the same as that of the Gross-Neveu model. All the modifications involve the couplings of the 
functions of p a ,o((3) to the other pi, j. After using the identities in the appendix, one finds that the 



net effect is to remove couplings between e aj o to e ft) i in the Gross-Neveu TBA (39), and replace 



them with a coupling between e a ^ to e Qi 2- The SU(N)/SO(N) TBA equations are 



N-l 



e . (g) _ T V/ (iV) !°° ^ - lnfl + e 6 ^^ 



=i 



h 31 27r 2cosh (iv(/? - 00/2) v + 



(43) 



The asymptotic conditions are the same as for the Gross-Neveu model. In fact, the only difference 
is that the second incidence matrix /(°°) is replaced with Z^°°\ These equations are displayed 
schematically in figure 4. 

Both cases can be conveniently summarized in the language of Dynkin diagrams: the Gross- 
Neveu model in figure 3 is described by (SU (N) , SU (oo)) , while the incidence diagram in figure 
4 for the SU(N)/SO(N) sigma model is described by (SU(N), SO (oo)). The latter TBA system 



was previously discussed in [35|, but without the association with the sigma model. As with 
all previous cases, one can check that the UV and IR limits of the TBA equations agree with 
known results, namely the central charge cuy = (N + 2)(N — l)/2 and the particles being in the 
representations 2p a . This computation in particular checks that these are all the particles in the 
spectrum, because additional (or fewer) particles would change this central charge. 

3.4 0(2P) Gross-Neveu models 

As with the models with SU(N) symmetry, I will start with the 0(2P) Gross-Neveu models 



[36], |37], |13[ (these are in fact the models Gross and Neveu originally studied). Like the SU(N) 
case, there are particles in all the fundamental representations with highest weights p a . This 
includes the spinor representations, which physically correspond to kinks. The mass spectrum is 
given by 

m a = m sin(a7r / (2P — 2)) mp-i = rnp 



2sin(7r/(2P-2)) 
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Figure 4: The incidence diagram for the SU(N)/SO(N) sigma model. There are N — 1 rows and 
an infinite number of columns. 

where the latter two correspond to the spinor representations. As opposed to the SU(N) case, 
for P > 4 there can be more than one representation with a given mass, as explained in detail in 
p3| . For any value of a there are particles in the representation [i a , but there may be additional 
ones as well. For example, for P = 4, there are particles in the vector and spinor representations 
(all 3 of them being 8-dimensional) of mass mi, particles in the antisymmetric representation 
(28-dimensional, weight /^) with mass m>2 = V / 3m-i, and a particle in the singlet representation, 
with mass mi- This apparently is related to representation properties of the Yangian; it turns 
out that the Yangian associated with 50(8) has a 29-dimensional representation, but not a 28- 
dimensional one. Under the 50(8) subalgebra of the Yangian, the 29 decomposes into 28 + 1. 
In the TBA equations below, the index a indicates all particles of mass m a , which presumably 



corresponds to an irreducible representation of the Yangian [15]. 

Luckily, the Bethe equations for 50(2P)-type systems were also found in [23, 24]. These were 
more or less conjectured based on analogy with the SU(N) case, but were proven up to some 
technical assumptions in |33]]. Basically, they amount to doing the computation by replacing the 
SU(N) incidence matrix J^*) with the 50(2P) incidence matrix I^ p \ The details for proving 
this are given in the appendix. The TBA equations for the 0(2P) Gross- Neveu models are 

p 



e -(B) - TTl {P) — lnfl + e^( 

aAI ) ~ h Loo 2vr cosh[(p - 1)03 - n v + 

X ^ l oo 2vrcosh[(P-l)(/3-/30] 1 H 1 + e J ^ 

These equations are displayed schematically in figure 5; the indices a and b now run over the 
nodes of a 50 (2P) Dynkin diagram. The correct central charge c\jy = P is obtained in the UV 



limit. This system was also discussed in [35]. 
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Figure 5: The incidence diagram for the 0(2P) Gross-Neveu model. There are P rows and an 
infinite number of columns. 



3.5 0(2P)/0(P) x 0{P) sigma models 

In || the 0(2P)/0(P) x 0(P) sigma models were shown to resemble the SU(N)/SO(N) sigma 
models discussed above. This is not terribly surprising, since the Lagrangian formulation of both 
is in terms of symmetric matrix fields. I will show here how their TBA systems are also similar. 

In H the exact spectrum and the S matrix S ss for the 0(2P)/0(P) x O(P) sigma models 
are found. Like the SU(N)/SO(N) case, there are particles are in all representations with 
highest weight 2/i a , where here a = 1 . . . P, although because of some peculiarities of the 0(2P) 
S matrices (and because of Yangian representation properties), there must be particles in some 
of the fundamental representations as well. The sigma model mass spectrum is the same as the 
0(2P) Gross-Neveu model, although of course the multiplicities differ. The TBA system for the 
0(2P)/0(P) x O(P) sigma models should not come as any surprise at this point. It follows from 
the 0(2P) Gross-Neveu model calculation just as the SU(N)/SO(N) calculation follows from 
that of the SU(N) Gross-Neveu model §. The TBA equations for 0(2P)/0{P) x 0{P) sigma 
models are 

aAP) - t[ ab J-oo 2n cosh[(P - 1)09 - /?')] \ + ) 

(oo) f°° d£ P-l / eaA p>)\ (45) 

h 2^ cosh[(P -!)(/? -/?')] \ +e ) [ b) 

The kernels and identities for this derivation are discussed in the Appendix. 



4 Massless sigma models with 6 = tt 

The results of the last section further confirmed the results of || for the S matrices of the 
SU(N)/SO(N) and 0(2P)/0(P) x O(P) sigma models when the instanton coupling 9 = 0. In 
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this section, I find the TBA equations for these sigma models when 9 = tt, further confirming 
results of 

The particles of the sigma models are massless when 6 = it. The reason is that both sets of 
models have stable infrared fixed points, the SU(N)i and 0{2P)\ WZW models, respectively. 
The S matrices for these flows were found in ||. Since the particles are massless, they are 
either left- or right-moving. Rapidity variables are still useful for parameterizing the energy and 
momentum of massless particles: E = p = me^ for a right mover, and E = —p = me~P for a left 
mover. The parameter m here is not the mass of the particle, but rather is the scale (analogous 
to Aqcd) which parameterizes the interactions. In condensed- matter language, it is the crossover 
scale. With these definitions, the rapidity difference is still an invariant in a collision. In a collision 
between a right mover and a left mover, the invariant is (E\ + E2) 2 — (pi +P2) 2 = m 2 e^ 1 ~^ 2 . In 
"collisions" between two right movers, the invariant is E1/E2 = e^ 1- ^ 2 . I put collisions in quotes 
because the S matrix is properly interpreted here as a matching condition on the wavefunction, 



as in fllPl). For more details on the S matrix approach to massless theories, see 0, 3C]. 

The spectrum and S matrices of these sigma models at 9 = tt are closely related to that of 
the corresponding Gross-Neveu model. For the SU(N)/SO(N) sigma model ||, 



s&M = S&09)/*2°G9) 
where X^ comes from fusing X as defined in (|33|): 

a b 

m = n n x 9 + [* + 3 - 1 - (° + h )/ 2 y N ) ■ ^ 
t=ij=i 

For N = 2, this reduces to the result of Qj. The reason for dividing out by X^\(3) in Slr is 
simple. Poles in Slr in the region < Im(/3) < tt are forbidden [Q], and all are contained in this 
factor. For similar reasons, the S matrices for the 0(2P)/0(P) x 0(P) sigma model at 9 = tt 
are || 

SfM = Si b L ((3) = S a G b N (P) 

?ab r o\ _ cab r o\ _ cab t o\ i y{P) , 



s&W = sf R (p) = s a G b N (p)/x^'(p) 



where 



XSF W = X aT~ 2) (m ( al P - 2) (** - 0) (47) 

and Sq N is the S matrix of the 0(2P) Gross-Neveu model. 

The TBA systems follow from the results in the last section, given the close relation with 
the Gross-Neveu models. The pseudoparticles are identical, so the densities p a j are labeled by 
two indices as before. However, in scattering, left movers stay left moving, and right movers stay 
right moving. Thus instead of densities p a ,o, now there are both p a< L and p a ,R- For the SU(N) 
case, the first of the Bethe equations (35) is replaced with the two equations 

JV-l N-l 
2TTP a M = ma e^+Y,Y^ ) *PbAP) + Y,( Y ab N) -^bS(P)+A^)*p b , L (P) 

6=1 6=1 

00 

^^ (oo) *p aj (/?) (48) 
i=i 
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AT— 1 N-l 



6=1 6=1 

oo 

r (oo) 



E a i ] *PM- (49) 
i=i 



The Bethe equations for the pseudoparticles (|36|) become 

iV-l oo 

27r Paj (/3) = af° ) * (p aji (/3) + Pa , R (P)) EE 4^ * * PM(^) (50) 



6=1 i=l 

Using the identities in the appendix gives the TBA equations 
N-l 

e a,j 



W = Yl^r^ * Infl + e^') 

^ ab 2vr 2 cosh(iV(/3 - /3')/2) I + 



=i 



V J*** r dP ' N In fl + e-^^l (51) 

^ ^ 7_ 2vr2coshfiVf/3-/3')/2) ln l 1 + e J (51j 



l=L,R,l...oo 



where j takes the values L, R, 1 . . . oo. These equations for the SU(N)/SO(N) sigma model at 
8 = 7r are identical to those for the SU(N)/ SO(N) sigma model at 8 = (f43|), once the labels are 
redefined (there j is takes the values 0, 1 . . . oo). However, that does not mean the solutions are 
the same. Because the 8 = theory is massive and the 8 = it theory is massless, the asymptotic 
conditions are different. Namely, as (5 — > ±oo, for the massive theory: 

eoo(/3 -> oo) — ► m a cosh(/3) 

while for the massless theory as /3 — > +oo 

tahifi oo) — > m a eP 

e a ji((3 — ► oo) — ► constant 

and as (3 — > — oo 

t-aLifi — ► -co) — ► constant 

£ai?(/3 -> -oo) — > m a e _/3 

The free energy ( |l~6| ) is modified in the massless case to 

FW(m,T) = -TE^a ^ [ e ^ln(l + e - £ ^)/ r ) +e -^ln(l + e-^^/ T )] (52) 

The equations for the massless theory are pictorially depicted in figure 6. 

The different asymptotic conditions do not affect the free energy in the ultraviolet limit 
m/T — * 0. Thus the free energy is the same in massive and massless cases, corresponding to 
that of a conformal field theory of central charge cjjy = (N + 2)(N — l)/2. This of course is 
the dimension of the manifold SU(N)/SO(N). In fact, because the TBA systems are identical 
except for the asymptotic conditions, the entire UV perturbation theory is identical in both 
cases. This is as it must be: instantons are a non-perturbative effect, and so the effect of the 



22 



I I 





Figure 6: The incidence diagram for the SU(N)/SO(N) sigma model with 9 = tt. There are 
N — 1 rows and an infinite number of columns. The cross-hatched circles represent eaL and e a R. 

instanton coupling cannot be seen in perturbation theory. Unfortunately, it is not known how to 
compute the perturbative expansion at non-zero temperature, except for the leading logarithmic 
correction The perturbative expansion at zero temperature can be computed explicitly by 
using a generalized Wiener-Hopf technique. This computation was done for the case at hand in 
H, and does indeed give the same results at 9 = and tt. 

On the other hand, the physics for 9 = tt is radically different from that at 9 = in the 
low-energy limit m/T — > oo. In the massive case the free energy in this limit is merely that of a 
dilute gas of massive particles, as in (|l7|) . However, the particles are massless when 9 = tt because 
the system flows to a non-trivial field theory in the low-energy limit. This flow is immediately 
apparent from the S matrix point of view, because the two-particle Lorentz invariant for a left 
and a right mover is oc in 2 , so the S matrix goes to a /3-independent constant value as m — > oo. 
The right-right and left-left matrices remain non-trivial, however, since the Lorentz invariant here 
is independent of m. Thus in the low-energy limit, the left and right sectors decouple from each 
other, but remain non-trivial. This is the behavior of a conformal field theory. The free energy 



must obey a relation like that of the UV limit, namely [20] 



vrT 2 

lim F= — c IR (53) 

m a — >oo b 

Here this gives cm = N — 1. This is the central charge of SU (N)x, confirming the flow discussed 
in§. 

In fact, since the left and right movers decouple in the IR limit, the TBA system for the 
right movers in this limit is obtained merely by removing the terms involving e a L from the 
equations. The resulting system is identical to that of the SU(N) Gross- Neveu model (|39|); only 
the asymptotic condition changes from e a o(/3 — > oo) — ► m a cosh/3 to e a ji(P — ► oo) — > m a e^. The 
TBA system for the left movers is the same, with the replacement (5 — > —(3. This close relation 
is a consequence of the fact discussed in B, that the effective field theory for the SU(N)/SO(N) 
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sigma model at 6 = tt in the low-energy limit is that of the SU(N) Gross-Neveu model at 
negative coupling. The sign change changes the sign of the beta function, meaning that while 
the Gross-Neveu model is an asymptotically-free massive theory, the critical point in the sigma 
model is stable. In another language, the different signs correspond to marginally-relevant and 
marginally-irrelevant perturbations respectively. 

Not surprisingly, the 0{2P)/0{P) x O(P) sigma model behaves in the same fashion. The 



TBA system in (45) applies to both massive and massless cases. Only the asymptotic conditions 
differ, as with the SU(N)/SO(N) model. As a consequence, the same cuv = P 2 is obtained for 
both 9 = and 9 = ir. In the massless case, the flow is to a conformal field theory with cjr = P, 
and the equations in the IR limit are those of the 0{2P) Gross-Neveu model. Thus indeed the 
flow is to the 0(2P)\ conformal field theory, confirming the results of ||. 

5 Perturbed coset models 

In |27], |38| it was shown how a G/H sigma model is related to a G^/Hi coset conformal field 
theory perturbed by a certain operator. In this section, I review this construction, and apply it to 
G/H = SU{N)/SO{N) and 0{2P)/0{P) x O(P). I find the exact free energy of the perturbed 
coset models, this approach shows promise for understanding whether other sigma models are 
integrable, as I will discuss in the conclusion. 

5.1 Perturbed coset models and sigma models 

A Gt WZW model is a conformal field theory with an infinite-dimensional symmetry algebra 
[39, 40]. This symmetry is an extension of a ordinary Lie algebra symmetry G. The symmetry 



currents are denoted J A (z) and J A (z), where A runs from 1 . . . dim(G). These currents satisfy 
the operator product 

J A (z)J B (w) = T ^ T2 + f ABCj °^ + . . . (54) 
[z — w) A z — w 

where the f ABC are the structure constants of the ordinary Lie algebra for G. The algebra (|54| ) 
is known as an affine Lie algebra or a Kac-Moody algebra Gk- The level A; is a positive integer 
for a compact Lie group G. The central charge (coefficient of the conformal anomaly) of the Gk 
WZW model is 

k dim G , . 

c = TT7T (55) 

where h is called the dual Coxeter number. It can be defined by fACDjBCD = h s AB / 2 . For 
G = SU(N), h = N, while for G = SO(2P), h = 2P - 2 (for P > 2). The primary fields of 
the WZW model correspond to representations of Gk- It is shown in |^0| that they have scaling 
dimensions 

^ = (Sy (56) 

where Cj is the quadratic Casimir defined by T A T A = Cjl, with the T A the generators of the 

Lie algebra of G in the jth representation and / the identity matrix. All the other scaling fields 

tAi 



arise from the operator product of the J (z) with the primary fields; it follows from ( p4[ ) that J 
has dimension 1 and therefore all fields have dimensions Xj plus an integer. 
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A coset conformal field theory Gk/Hi is formed from a Gk WZW theory and a subalgebra 
Hi. The energy-momentum tensor is constructed from the generators of Gk not in Hi [41]. The 



central charge of this new conformal field theory is c(Gk) — c(H{). The level / of the subalgebra 
H is determined is given by I = kr, where r is a group-theory factor called the index of the 
embedding of G into Hi. For the embedding of SO(N) into SU(N), r = 2 (r = 4 for N=3), 
while for the embedding of O(N) x O(N) into 0(2N), r = 1 (r = 2 for N = 3). 

The fields of the G^/Hi conformal field theory are constructed by decomposing a field 4>q in 
Gk into representations of the Hi subalgebra. Because the energy-momentum tensor obeys the 
orthogonal decomposition Tq = Th + T g / h , the decomposition of <j>G must be of the form 

<f>G = ®a(f>G/H®<f>H- ( 57 ) 

The coefficients (j> G / H of this decomposition are the fields of the coset model Gk/H\. 

These coset conformal field theories a priori have nothing to do with G/H sigma models. 
The former are massless, and do not have a global symmetry G, while the latter are gapped with 
a G global symmetry. Thus for the two to correspond, the coset model must be perturbed by 
some operator. Moreover, the coset model has a G global symmetry when k — > oo. These and 
other considerations led to a conjecture made in p8| . This conjecture is that the sigma model for 
G/H is equivalent to the k — ► oo limit of the coset conformal field theory perturbed by a certain 
operator. The operator is obtained by using ( |57| ) to decompose the currents J A into fields in 
Gk/H[. For the cases of interest here, G/H is a symmetric space, meaning that there is no normal 
subgroup of G containing H other than G itself. A consequence of G/H being a symmetric space 
is that the generators of G not in H form a real irreducible representation of H ]^2|] . Thus when 
a field J A (z) is decomposed into representations of H in ( |57| ) there is only one term on the right 
hand-side. The resulting field in Gk/Hi is denoted by J A . The fields J A form a real irreducible 
representation of H, of dimension cjjv = dimG — dimH. The operator O a is defined as 

o = j2j a w a ^)- ( 58 ) 

A=l 



The conjecture of [38] can now be stated precisely: the G/H sigma model is equivalent to the 
Gk/Hi coset conformal field theory perturbed by the operator O a in the limit k — * oo. 

The conjecture passes a few simple checks. The ultraviolet limit is obtained by removing the 
perturbation of the coset model. From ( |55| ) it follows that the central charge of the Gk/Hi theory 
as k —* oo is indeed c\jy = dim G— dim H as in the sigma model. Moreover, when one decomposes 
J A into representations of H for A in G but not H , the resulting field <pjj has dimension going to 
zero as k — > oo, because the quadratic Casimir in ( f)6|) is independent of k. Thus the field J A has 
dimension 1 in this limit, so the perturbation O a is of dimension 2 and so is naively marginal. It 
is not exactly marginal - this is the phenomenon of dimensional transmutation and asymptotic 
freedom. Therefore the coset and its perturbation have the general properties of a sigma model. 



Further support for this conjecture is discussed in [38]. For example, it has been shown to be 
true for the principal chiral models |^3|], and in the sphere sigma model p7|| . The results in this 
section give strong further evidence in support. 

The models of interest in this paper are the SU(N)k/ SO(N)2k and the 0(2P)k/0(P)k X 
0(P)k conformal field theories perturbed by O a . The former theories have 

k(k-l)(N + 2)(N-l) 
CUV ^ N) = (N + k)(N-2 + 2k) 
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while the latter have 

c uv {k,P) - 



(P -2 + k)(2P -2 + k) 

To find the dimensions of the perturbing operators requires a little group theory. Fields in the 
adjoint representation of SU(N) decompose under the SO(N) subgroup as 

The representation of dimension N(N — l)/2 consists of the generators J A with A in the SO(N) 
subgroup as well. Thus the operators J A are in the symmetric representation of SO(N)2 k , of 
dimension N(N + l)/2 — 1. The quadratic Casimir of this representation is C syrn = N. Since 
the dimension of J A is always 1, the dimension of O a in the SU{N) k / SO{N)2 k conformal field 
theory for TV > 2 is 

2N i k-1 

Xa ~ 2 ~ N - 2 + 2k ~ N - 2 + 2k 
Similarly, the adjoint representation of 0(2P) decomposes into 

( p {2F _ 1)) ^^, 1 Wi,^U(p,p) 



under the O(P) x O(P) subgroup. Thus the operator O a here is in the (P,P) representation of 
0(P)k x 0(P) k . The quadratic Casimir of the vector representation of O(P) is (P — l)/2, so 

P- 1 k-1 

As far as I known, these perturbed conformal field theories have never been studied in the 
literature. 



The role the instanton coupling 9 takes in the conjecture of [27, [38| is quite interesting. The 



action of the perturbed conformal field theories can be denoted schematically as 

S = Scft + X J d 2 zO{z,z). 

It follows from simple scaling considerations that mass scale m in the theory is related to A 
by m oc lAI 1 ^ 2-2 ^. If the theory has a Z2 symmetry under which O — > — O then the theories 
with positive and negative A are identical. In general, they are not. A well known example is 
the SU(2)k x SU(2)i/SU(2)k+i "minimal" models of conformal field theory perturbed by O a 
(usually called (^13 in this context). With one sign of A, the model is massive. With the other 
sign, the model flows to the minimal model with k — 1 f4^| , so the excitations are massless. In the 
SU(N) k /SO(N) 2k and 0(2P) k /0(P) k x 0(P) k cases for k > 2, the two signs of A give different 
theories as well, one massive and the other massless. In the k — * 00 limit these differing theories 



correspond to 6 = and 9 = ir respectively. This was argued in 27] for SU{2) k /0{2). Strikingly, 
one can also see from the perturbed conformal field theories here that that the different sign 
affects perturbation theory only at the order X k . Thus as k — * 00, the different sign does not 
affect perturbation theory. Its only effects are non-perturbative, just as they must be if the 
change A — > —A is to describe the effects of a 9 term. 
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5.2 The particle spectrum 

Here I discuss the particle spectrum of the perturbed conformal field theories just defined. 

The results for the simplest cases k = 2 are already well known (when k = 1, the models are 
trivial). The SU(N)2/0(N)4 : conformal field theories are known as the Zn parafermion theories; 
the equivalence to the better known coset description SU(N)\ x SU{N)\/ SU{N)2 was shown in 
p5| . The perturbation O a of dimension 2/(N + 2) is called the thermal operator here. This is an 
integrable field theory, with S matrices derived in The spectrum consists of iV — 1 particles, 
with mass [^] 

m a = man 

This is the same mass spectrum as in the SU(N) Gross-Neveu model and the SU(N)/ SO(N) 
sigma models discussed above. The degeneracies are different: there is only one particle of each 
mass in the parafermion model, while in the other cases, there are multiplets of particles in 
SU(N) representations with highest weights \x a and 2[i a respectively. The parafermion theory 
has a Zjy symmetry, but no SU(N) symmetry. 

Likewise, the 0(2P)2/0(P)2 X 0{P)2 conformal field theories are the T>2p parafermion theo- 
ries. Their symmetry group is not 0(2P), but instead the dihedral group T>2p- The equivalence 
to the usual formulation 0{2P)\ x 0{2P)\/0{2P)2 formulation of these parafermion theories 
can be shown using the techniques of p5] . These theories have c = 1 for any P. The pertur- 
bation is of dimension 1/P, and so the massive theory corresponds to the sine-Gordon model 
at f3 2 = 8vr/P. This is of course integrable 111]], and in fact corresponds to the "reflectionless" 



points of sine-Gordon, where the scattering is diagonal. The spectrum consists of P particles, of 
masses 

Tfl 

m a = msin(a^/(2P - 2)) m P ^ = m P = 2 sin( ^ /(2p _ 2)) 

The particles of masses P and P — 1 are the kink and antikink of the sine-Gordon model. This 
mass spectrum is the same as that in the 0{2P) Gross-Neveu model, and the 0(2P) /0{P) x 0{P) 
sigma model, but with multiplicity 1 here. 

The fact that the mass spectrum of the k = 2 perturbed coset models are the same as the 
corresponding sigma models is already a strong piece of evidence in support of the conjecture 
of l38|]. The issue now is to find the spectrum and S matrices for general k. For N = 2 and 



P = 2, the answers are given in [48], but otherwise these models have not been discussed in the 
literature. I will solve this problem for all k. 

To understand the particle spectrum in an integrable model, it is crucial to understand the 
symmetries of the model. For the sigma models, this symmetry algebra is an ordinary Lie algebra 
G. I conjecture that the perturbed coset models are invariant under a one-parameter deformation 
of G called the quantum- group algebra U q {G). The particles in the perturbed G^/Hi models form 
finite-dimensional representations of U q (G), with the parameter q = e l7T /( k + h ) ; where h is still the 
dual Coxeter number. As k — > oo, q — > 1 and the algebra reverts to the usual G Lie algebra. 
All known integrable perturbations of coset conformal field theories are proven or believed to be 
invariant under some such quantum-group algebra. For example, for models where G = H x H, 



this was discussed in detail in [43|. For other models, this was discussed in [49]. For the cases of 



interest here, the particles in the SU (2)^/0(2) were shown to form a representation of U q (SU (2)) 



m 



To give a concrete example, U q (SU(2)) is the algebra 



[S z , 5±] = ±2S±, [S+, S-] = ^ q J Sz (60) 
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When the parameter q = 1, this reverts to the usual SU (2) Lie algebra. A nice physical realization 
of this algebra is discussed in [50], where it is shown how when the Heisenberg spin chain is 
deformed into the XXZ spin chain, the SU(2) symmetry is deformed into U q (SU(2)). The 
properties of the representations of U q {G) can be quite different from those of G when q is a 
root of unity other than 1 or —1. For example, the right-hand-side of the last equation in ( |60| ) 
vanishes on states with 2S Z = p when q v = 1. This means that representations with maximum 
value of 2S Z greater than p are reducible. In other words, the only irreducible representations 
have \2S Z \ < p, as opposed to ordinary SU(2), where there are irreducible representations with 
any integer value of 2S Z . 

Particles in a representation of a quantum-group algebra are most conveniently treated as 
restricted kinks [|5l] l. Consider a field 4>, with a potential V((f>) tuned so that there are degenerate 
minima, which I will sometimes call vacua. Then kinks are field configurations with 4>(x = 
— oo) one minimum of the potential, <j)(x = oo) another. The kinks in the perturbed coset 
models form what are called "restricted solid-on-solid" , or RSOS, representations of the quantum- 
group algebra. The name comes from the statistical mechanical lattice models in which these 
representations first arose [52]. 

For U q (SU(2)), these restricted kinks are easy to describe. They interpolate between the 
minimum of a potential which has k + 1 minima in a row. For example, the potential V{<f>) = 
(j) 2 ((f> 2 — l) 2 has three minima at (f> = 0, ±1; the potential V{4>) = (4> 2 — l) 2 (</> 2 — 9) 2 has four 
vacua in a row. Kinks in these sorts of potentials provide representations of the quantum-group 
algebra U q (SU(2)) with q a root of unity. The two-dimensional representations are kinks which 
interpolate between adjacent vacua. Such representations behave just like ordinary SU(2) spin- 
1/2 representations. For example, for k = 2, there are three minima labeled 0, ±1, and the 
generators S± exchange the states </>(— oo) = and </>(+oo) = ±1. To construct the larger 
representations, one can take the tensor product of smaller representations. The rules are just 
like that of ordinary SU(2): for example, the tensor product of two spin-1/2 representations 
decomposes into the sum of a spin-1 and a spin-0 representation. The one catch is that for q a 
root of unity, the larger representations are reducible. For k = 2, spin 0,1/2 and 1 are all the 
irreducible representations. This is clearly apparent from the kink picture, because for k = 2 
there are only three vacua: an irreducible spin-3/2 representation requires four vacua. Moreover, 
even the allowed kinks are restricted. Restricted means that multi-kink configurations must obey 
the rules implied by the potential. The number of AA-kink states is much less than the number 
of one-kink states to the Af^ 1 power. For k = 2, in fact, there is only one way to construct a 
multi-particle state from spin-1 particles: the vacua must alternate between +1 and —1. The 
restriction is so strong that the kink structure gives no new degrees of freedom, so it can be 
viewed as a normal particle. 

Perturbed coset models with restricted kinks are already widely known. The SU(2)f : x 
SU (2)i I SU (2) k+\ minimal models perturbed by O a are integrable. The particles are spin-1/2 
U q {SU{2)) kinks, where q = e in ^ k+2) ||, g|. There 

are thus k + 1 vacua here, with the kinks 
interpolating between adjacent vacua. The k = 1 case corresponds to the thermal perturbation 
of the Ising model (free Majorana fermions). Since there are only two wells when k = 1, all the 
kink can do is go back and forth, and one can forget it is a kink. For the case SU(2)i t /0(2), the 
particles are spin-1 U q (SU(2)) kinks The k = 2 case here also corresponds to the thermal 
perturbation of the Ising model. In this description, there are three vacua, but the kinks are of 
spin 1, so again all they can do is go back and forth: there is only one state for a given number 
of particles. 

For general algebras U q (G), the restricted-kink structure is more complicated. The potential is 
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defined so that the minima correspond to the highest- weight states of the quantum-group algebra 
allowed at that value of q. For simply-laced algebras, the allowed weights Yl a °a^a must satisfy 
^ a c a < k. In this language, for U q {SU{2)) with q 4 = 1 (k = 2), the three minima correspond 
to highest weights 0, /J,i,2fj,\, where fj,% is the sole fundamental weight of SU{2). The kinks form 
representations of the algebra, so each kink is also labelled by a weight. The rule is then that 
there can be a kink of representation r a interpolating from the vacuum 7 to the vacuum 5 if the 
corresponding representations obey the tensor product 



r a <g> r 7 = r<5 © . . . 

I said "can be" because it depends on the specifics of a given theory if such a kink actually does 
appear in the spectrum. There are a number of subtleties with this picture for general groups 
and representations, but it is not necessary to understand them for this work. 

Given a particle spectrum consisting of restricted kinks, the S matrix can be found using the 
Boltzmann weights of the corresponding lattice statistical-mechanical model, which is usually 
known as the R matrix. For models with particles in the fundamental representations, this was 



discussed in [51, 43, |53[. I emphasize that by corresponding lattice model, I do not mean a lattice 
model whose continuum limit is described by a field theory with this S matrix. I mean that there 
is some integrable lattice model whose Boltzmann weights are proportional to the S matrix. In 
the corresponding lattice models, the variables which placed on sites of the lattice play the role 
of the vacua, while the kinks correspond to the states on the links. The rapidity difference in the 
S matrix corresponds to the spectral parameter in the lattice model. The scattering of kinks in 
representation a from one in representation b is given by the matrix S ab oc R ab (as before, a and b 
are not the matrix indices, but rather label the different matrices). The prefactor is not of interest 
to the lattice model, since it merely multiplies the partition function by an overall factor. It is 
of course of great importance to the S matrix theory. The R matrices for the RSOS models are 
trigonometric solutions of the Yang-Baxter equation. They can be written in the form (|5|), where 
the f£ b are trigonometric functions (as opposed to the rational functions appearing in the sigma 
models). They are given explicitly for the fundamental representations of all the quantum-group 
algebras U q (G) in generalizing the SU(2) results of p^| . The fusion procedure also can be 



used to construct the R matrices for kinks in the representations 2/j, a [p5 |. 

The spectrum of the perturbed coset models is easy to obtain, given the sigma model result. 
The kinks must be in the same representation of U q (G) as the particles are of G. For example, 
for the case SU(2)k/U(l), the kinks are in the spin-1 representation of U q (SU(2)), while the 
particles in the sigma model are in the spin-1 representation of SU{2). When k is finite, the 
vacua are restricted, but the restriction is removed as k — > 00: the particles in the sigma model 
no longer need be viewed as kinks. Similarly, for the massive perturbed SU{N)k/ SO{N)2k and 
0{2P)k/0{P)k x 0(P)k models, the kinks are in all representations 2[x a for a = 1 . . . N — 1 and 
a = 1 . . . P respectively. The vacua are all weights ^ a c a fi a with ^ a c a < k. 

Note also that the H Gross-Neveu model is obtained by taking k — ► 00 in the perturbed coset 
models H k x H x /H k+l M, || M. 



5.3 The free energy of the perturbed coset models 

The derivations of the TBA equations for the perturbed coset models requires diagonalizing the 
transfer matrices formed from the kink S matrices. The computation is very similar for those of 
the sigma models, because the analysis of [23, pi] applies to the RSOS models. 
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It is simplest to first discuss the case k = 2, where the perturbed coset models reduce to the 
well-studied parafermion theories. As explained above, the kink structure is trivial: there is only 
one particle for each representation a = 1 . . . N— 1 or a = 1 . . . P. The scattering here is diagonal 
but non-trivial. The S matrix element for scattering a particle of type a from one of type b for 
SU (N) parafermions is 

S ah (P)=X%\(3) 

where xj£\(3) is defined in (f46|). For the 0(2P) parafermions, the S matrix elements are 

(P) 

^ab \P)i as defined in (|47| ). The TBA equations instantly follow from using these S matrices 
to give the kernels in fllq). There are no pseudoparticles because the scattering is diagonal, so 
the only functions which appear can be labelled e a ^. Note the distinction with the sphere sigma 
model, where the only functions which appear are e\j in the present notation. For the SU(N) 
parafermions, the TBA equations are [p7f| 



N 1 /'OO r l /Df 

e a ,o(P) = m a cosh - Tin (l + e " e - °W' T ) + T £ / ^-A^\(3 - ff) In (l + e 

b=i J -°° Z7T 



where A^P is the same kernel which appeared in the Bethe equations above, and is given explicitly 
in the appendix. This can be simplified greatly by using the fact that A and K are inverses, 
giving 

e aM -1 2^1 ab J ^ ^ cosh[N{(3 _ 0/)/2] ln^l + e ) (61) 

where the asymptotic condition e aj o — * m a cosh (3 as (3 — > oo is implied. The incidence matrix 
couples only "adjacent" functions; it is displayed by restricting the diagram in figure 3 or 4 to 
have only one column. For the 0(2P) parafermions, the kernel A^ is replaced by A.^ [47]. 
This results in the TBA equations 

Thus the TBA equations for the k = 2 cases amount to those of the corresponding sigma models 
with all the pseudoparticles removed. 

The TBA equations for general k are also found by truncating the equations for the corre- 
sponding sigma model. The reason is simple to describe schematically. Each irreducible repre- 
sentation of the quantum-group algebra is associated with some transfer matrix. Relations like 
the fusion relation (^) relate the different transfer matrices. The fact that there are only a 
finite number of irreducible representations of the quantum-group algebra means that the fusion 



relations relating all these transfer matrices truncate 23, 24]. In the Bethe ansatz equations, 
this means that there are only a finite number of pseudoparticles. In the TBA equation, the 
index j in the functions e a j now runs only from . . . k — 1 in the massive case. 



This derivation of the Bethe equations is covered in detail in 24, 34]. For the SU(N) 
case, for example, the Bethe equations for the pseudoparticles are very similar to (|36|), but are 
modified to 

JV-i fe-i 

2^ Pa M = of * p a M - Y, A f * ^ * ( 63 ) 

6=1 1=1 
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The equations for the 0{2P) case are modified in a similar fashion. The S matrix prefactor is 
modified as well; the kernel for the SU{N) k / SO{N) 2k case is given in the appendix. 

The result of these modifications is that the TBA equations are truncated. Like the TBA 



equations (39,43,i~i|]4"5|) they are of the form 



rank G r s 
T ]T Q ab 



6=1 



h 



k-1 

1=0 



2vr 2cosh[/i(/3-^')/2] 
00 dp h 



2vr 2cosh[/i(/3 - f3')/2 



In ( 1 + e' 6 *^ 



(64) 



where h is the dual Coxeter number for G, which is N for SU(N), and 2P — 2 for 0(2P). The 
rank of SU(N) is N — 1, and the rank of 0(2P) is P. The matrices Q and R are all incidence 
matrices. For the various models considered here, the results are given in the following table. 



coset model perturbed by O a 


Q 


R 


behavior when k — > 00 


SU(N) k x SC/(iV)i/S£/(iV) fc+1 


I(N) 


/(*) 


SU(N) Gross-Neveu model 


SU(N) k /SO(N) 2k 






SU(N)/SO(N) sigma model 


0(2P) k x 0(2P) 1 /0(2P) k+1 




l(k) 


0(2P) Gross-Neveu model 


0(2P) k /0(P) k x 0(P) k 




X {k) 


0(2P)/0(P) x 0(P) sigma model 



In all cases, the usual asymptotic conditions apply. All the TBA equations in this paper are 
contained in this table. One can check that the central charges resulting from taking the UV 
limit of the TBA equations are indeed those of the corresponding conformal field theories for any 
value of k. This is an enormous check on all the results of this paper. 



5.4 Flows between coset models 

I showed for the sigma models that the TBA equations for 8 = and 9 = it are identical, with the 
only difference being in the asymptotic conditions. The same behavior should happen for the two 
signs of A in the perturbed coset models (^). The TBA results for the perturbed coset models 
make it possible to understand the flow when the perturbation is massless. The TBA equations 
(|64| ) and the table still hold, except that the asymptotic conditions given in section 4 apply here. 
The sum over I now runs from L, R, 1 . . . k — 2. The IR fixed point can be read off from the 
equations, as described above for the sigma models. Removing say the left moving particles from 
the SU(N) k /SO(N) 2k perturbation gives the diagram for the SU{N) k -i X SU(N)i/ SU{N) k 
models. Thus the flow is between the conformal field theories 

SU(N) k SU(N) k ^ x ST^AQi 
SO(N) 2k ~^ SU(N) k 

Likewise there is a flow 

Q(2P) k Q(2P)k-i x Q(2P) l 

0{P) k x 0(P) k ~^ 0{2P) k 

As far as I know, these flows were previously unknown. When k = 2, there is no flow: the two 
cosets are already equivalent. By using the equivalences between different coset models derived 
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in [45], one can described these flows in different ways, if desired. For example, the latter also 
amounts to a flow 

Q{k) P x Q(k) P Q(2P) fc _ 1 x Q(2P) 1 
0(k) 2 p ' 0{2P) k ' 

Going backwards, one can read off the spectrum and particles for these massless perturbations. 
The kinks must be massless, and in all representations /U a , and are either left or right-moving. 
The vacua correspond to all weights ^ a c a/ u a with Y^, a c a < k — 1. This shift of k — > k — 1 indicates 
the quantum-group parameter q is different for the massless and massive perturbations, but I do 
not know the reason for this. In the coset models x Hi/H^i there are two quantum-group 
symmetries for both perturbations presumably the same thing happens here. 



6 Conclusion 

In this paper I have described how to compute the exact free energy in integrable two-dimensional 
sigma models. This definitively establishes that when 9 = tt, there are non-trivial fixed points 
for two sets of sigma models. It also yields the exact free energy and susceptibility when 9 = 
and when 9 = tt. 

The big open question is if other sigma models are integrable. The grail in particle physics 
is probably the CP N ^ = SU(N)/SU(N - 1) x 17(1) models. They have been widely studied 
because they allow instantons and are tractable in large N. (The models studied above have a 
parameter N and have instantons, but they are difficult to treat in large N. The reason is that 
they are matrix fields: the number of fields at large N grows as N 2 , not as N.) In particular, the 
CP N models allowed Witten to conclude that instantons were not important in real-world QCD 



[57]. It would be very interesting to prove Witten's results directly, instead of relying on large 
N. 

Virtually all the symmetric-space sigma models have arisen in various condensed-matter ap- 
plications [|5Sf| , but the grail here is the U(2N)/U(N) x U(N) "Grassmanian" model. The reason 
is that in the replica limit N — > 0, this is believed to describe the transition between quantum 
Hall plateaus j|] . This transition is experimentally realized, and good numerical and experimen- 
tal measurements have been made of critical exponents. These critical exponents should arise in 
some conformal field theory, but it is still not known which one. Solving the sigma model as a 
function of N would presumably solve this problem. 



So why are sigma models integrable? In some sigma models (see e.g. [59, |5f|), one can 
find non-local conserved currents. Although the existence of non-local currents does not prove 
integr ability, it is a good indicator. Often these non-local currents are often associated with 
quantum-group or Yangian symmetry algebras. In the 0(N)/0(N — 1) models, one can prove 



the non-local currents of [59] are the generators of an infinite-dimensional symmetry algebra 
called the Yangian This proves the integrability of these sigma models. Unfortunately, this 
result has not yet been extended to other sigma models. 

So are other sigma models integrable? An old result (see e.g. |56|]) suggested that the only 
integrable symmetric-space G/H sigma models are those where H is a simple Lie group. The 
reason is that they found that the non-local conserved currents coming from the classical sigma 
model (the limit of g small) are not conserved once loop corrections are included. This certainly 
does not prove the model is not integrable, because it is possible that some or all of the classical 
conserved currents can be modified so that they are conserved in the full theory. 

A simple Lie group has only one factor. Thus the symmetric spaces with H simple are 
0(N)/0(N-1), SU(2N)/Sp(2N) and SU(N)/SO(N), and the principal chiral models HxH/H. 
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All of these models are indeed integrable. However, the 0(2P)/0(P) x O(P) models are also 
integrable, but H is not simple! Thus the suggestion of [^f| is in not true here. It is not clear 
whether this is a fluke of this model, or other symmetric-space sigma models are integrable as 
well. It would be most interesting to construct the non-local conserved currents here explicitly, 
to understand how they remain conserved even in the full theory. Some interesting results for 
the classical model were found in but they await generalization to the quantum case. 

The S matrices described above are all what are known as rational solutions of the Yang- 
Baxter equation. This means the S matrices are rational functions of the rapidity (except for the 
prefactor). Yangians are all associated with rational solutions of the Yang-Baxter equation, so 
the results described above certainly imply that there is a Yangian symmetry in all the integrable 
sigma models. In fact, this is the reason for the extra particles in the models with 0(2P) 
symmetry. The representations of the Yangian of 0(2P) are larger than that of its subalgebra 
0(2P). The particles at a given mass are in a reducible representation of 0(2P), but in an 
irreducible representation of the Yangian. This poses an interesting question: is there any way of 
telling which representations of the Yangian yield the particles and S matrices for an integrable 
field theory? And if so, what are these theories? Unfortunately, the technology of Yangians does 
not seem developed enough yet to answer these questions. 

The results of [3£| discussed in section 5 do suggest an alternate approach to finding integrabil- 
ity in sigma models. It is much easier to look for conserved currents in perturbed conformal field 
theory than it is in sigma models. For example, it was noted in |3£| that there are (at least to low- 
est order in perturbation theory) conserved non-local currents in the SU(N)k/SU(N — 1)& x U (1) 
coset models perturbed by the operator O a . Thus one expect these currents to persist in the 
CP 1 sigma model, obtained by taking k — > oo. Even if these currents do remain in this limit, 
this does not prove the CP N_1 models are integrable. However, at the very least it would indicate 
that interesting behavior in the sigma models is still lying yonder. 



My work is supported by a DOE OJI Award, a Sloan Foundation Fellowship, and by NSF 
grant DMR-9802813. 



A Kernels and identities 

A.l SU{N) 

One set of kernels I use comes from the prefactors of the S matrices. These kernels are defined 

as 



These kernels arise in the prefactors of the SU(N) parafermion theories, the SU(N) Gross- Neveu 
models, the SU(N)/SO(N) sigma models, and the SU(N) S / SO(N)2 S perturbed coset models, 
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respectively. The reason for the extra factor in the definition of A^' will become apparent below. 

The kernel appearing in vector- vector scattering is defined as 7^ = Y±i • It is most useful to 
give the kernels in Fourier space. To make the equations look a little nicer, I define the Fourier 
transform with normalization 

/>) = - r e^'-m (65) 

I use this definition of Fourier transformation for any kernel in a model with SU(N) symmetry. 
A fact useful for obtaining the TBA equations in this paper is that if f(ui) = l/cosh(u;) then 
f(/3) = jV/(4vrcosh(A^/?/2)). 

For the Gross-Neveu models, by using the S matrices in fn] , 32] one finds after some after 
some manipulation [34] 



r(N) f/ A _ x , _ >| siah((N ~ sinh(frq; 

r( N ) _ rp n j ,i W„„1 D jpab 



Kb ( w ) = °ab - e 1 1 . , ( , T >,.,/>, (66) 

ao sinh(JV(j) smh(a)J 



for a > b, with lj, = KV To find the kernels F appearing in the SU(N)/SO(N) sigma 
models requires even more work. Using the results of || for the S matrices, I find 

Z {N h >(*) = S ab - e -|H 4cosh H S mh((jy-aV)sinhM 
Notice how the Fourier transforms are related: 

r{N)r, A A . _ „-2u^to,,\ (vW,,.,> . )iih 



Z£>(u) - S ab = e~*> sinh(2u,) [Y£>{u) - 5, 

This relation is useful in proving various identities. Finally, for the perturbed coset models, one 
has 

(jv, s ), > 4cosh(w) sinh((iV — a)u) sinh(6w) sinh((s — 

Qb ' ( U ) = "ah ■ wju \ ■ . / \ 

smh(iVu;) smn(sw) 

Note that Z^ N ' = ^( Ar '°°) ) in accord with the idea in p7| , |38| that the sigma models can be 
obtained as the limit of perturbed coset models. 

The kernel A^f) arises in several places. The functions X^jp are the S matrix elements for 
the SU(N) parafermion theories, and appears as part of the prefactor in the Gross-Neveu and 

iW 



SU(N)/SO{N) sigma models. A$ also arises in the Bethe ansatz diagonalization. It is 



- (s) 2sinh((s - j)w) cosh(w) sinh(L;) 

A )l ( u ) = ■ u ( \ ■ Z ( ^ ( 68 ) 

for j > I, with Ay = A^ . Other kernels arising in the Bethe ansatz diagonalization are 



in the Gross-Neveu models, and 



8f ( W ) = (69) 



?{ ) = 2sinh(( S - J >)cosh(a,) _ 
J sinh(su;) 
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in the sigma models. Notice that r and a are related via ©. Naively, this seems to imply 
Tj S \lo) = 2 cosh.(u))dj S \uj) , but this is not quite true. The 5j± appears in (70) after a careful 
analysis of the Fourier transforms; note that the correct forms vanish as uj — > oo. 

(s) 

The inverses of the matrices are very useful. By using the Fourier transforms, it is simple 
to derive the identity 

y>,?*4?(/?)=w 



^-jk *^ki VP) = °\P)°ii 
fc=l 



where 



j{s) 

Kf(u) = 8 jl - - s ±- - (71) 
• Jt J 2cosh(w) 

where I-f is the incidence matrix for the algebra SU(s), defined in (|38|). More generally, the 



incidence matrix for a simply-laced Lie algebra is twice the identity minus the Cartan matrix, 
and is conveniently pictured by the Dynkin diagram. I denote the incidence matrix for SO (2s) 
as . Other useful identities are 



and 

s-l 



Useful identities involving the £ matrix prefactors are 



and 



Z^(u)-5 ab = A^(u) 



2 cosh(cj) 

.(oo). 



2 cosh(w) 

The extra Sji in (f7C|) is crucial to obtaining the right identities. 



A.2 0(2P) 

The S matrices and prefactors for the 0(2P) Gross-Neveu models are given in [jl3j, and those 
for the 0(2P)/0(P) x O(P) in §. The kernels are defined as and Z$ respectively. The 
Fourier transform used below is that of ( |65| ) with N replaced with 2P — 2. 

For a,b = 1 ... P — 2, the Gross-Neveu kernels are closely related to the SU(2P — 2) kernels, 
namely 

V^(,A - y(2P~2)( ,,\ , y-(2P-2) / A 
Sab \ U ) - Y ab (CO) + y 2P - 2 -a b\ U ) 

j. |^| cosh((P — 1 — a)uS) sinh(6u;) 

cosh((i- > — 1)uj) sinh(u;) 
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^( P) ^(1P 2) 

for a > 6, with y} = . For those involving the spinor representations s and s (the nodes 

labelled P — 1 and P), the kernels are 

- f {P) -V (P) - 1 >l smh(Po;) 



2sinh(2u;)cosh((P- l)w) 



v( p ) - - „M sinh(ao;) 



aP 'aP-i 2sinh(cj)cosh((P- l)w) 

where in the latter a = 1 ... P — 2. 



The result of [23, 24] for the Bethe equations for 0(2P — 2) says that the equation for the 
eigenvalue (|34T) and the first Bethe equation (55[) are 



oo 

2vrP a , (/3) = m a cosh /? + ]T ) * p 6>0 (/3) - E ^ {oo) * £,,,(0). (72) 



6=i ;=i 



These are virtually identical to those for SU(2P — 2), with Y^' replaced by Y^ P 2 \ In partic- 
ular, the kernel CT ( 2P - 2 ) is still given by (|69|). The other Bethe equations are now 



P oo 

2K Pa M = af* * PaM EE 4^ * 4? * P*,i09) (73) 



6=1 1=1 

where 



£W(a;) = Sj i- P -{ P r[ I (74) 
^ v 1 \ J 2cosh(u;) I 

(P) p— — j 

where J-, is the incidence matrix for the algebra SO(2P), defined above in (B2). The reason for 



j7 

the P — j and P — I indices is that above it was convenient above to define the spinor nodes as 
and 1, whereas here I have defined them as P and P — 1. 

The proof of the TBA equations is now basically identical to that done for the SU(N) Gross- 
Neveu model. The reason is that the kernels here satisfy basically the same identities as the 
SU (N) case. Namely, one can define the matrix inverse A of /C, just like A is the inverse of K. 
One finds that 

-i±\nX a \p)=5 ab 5{P)-A {P \p) 

Then 



2 cosh( 



CO 



Using this and the identities in the first appendix gives the 0(2P) Gross- Neveu TBA equations 
in ((Op. 

For the 0{2P)/0(P) x O(P) models, the proof is the same as for the SU{N)/SO{N) models. 
The only new identity needed is 

m f ,.,\ _ x . _ A( p )r,^ ( y 2° o) (^) _ ,\ 



2 cosh(w) 



From the prefactor given in it follows that this identity holds for a = b = 1. However, I 
have not been able to prove it in general. The reason is that the S matrices for particles in the 
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representations 2/j, s and are not known explicitly, so it has not been possible to work out 
the prefactors involving these particles. However, I have checked that if they obey the above 
identity, then they are consistent with the massive TBA. By consistent, I mean that the TBA 
equations are the same as in the massive case with only different asymptotic conditions, so that 
the perturbative expansion of the free energy is the same for 9 = and = ir. I have also checked 
this consistency for the energy at zero temperature in a magnetic field, extending the analysis of 
H to the particles in representations fi s for the massless case and 2fi s in the massive case. 



As a tangential note, the 0(2P)/0(2P — 1) sigma models are integrable as well [11]. Their 
spectrum consists of a single multiplet of 2P particles in the vector representation, with no bound 
states. The TBA equations are very similar |p^| , but a in p a Q can only be 1. The other e a j still 
have a = 1 ... P. Because there are no bound states, the prefactor J 711 is not the same is in the 



Gross- Neveu models: X{0) needs to be removed from the prefactor [11 1. The kernel appearing in 
the TBA equations is therefore yfPifi) - 5 ab 5{f3) + A^(J3). Using this with the ab ove Bethe 
equations gives the TBA equations given in p8| . 
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